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Abstract 

I consider a one dimensional system of particles which interact 
through a hard core of diameter a and can connect to each other if 
they are closer than a distance d. The mean cluster size increases as 
a function of the density p until it diverges at some critical density, 
the percolation threshold. This system can be mapped onto an off- 
lattice generalization of the Potts model which I have called the Potts 
fluid, and in this way, the mean cluster size, pair connectedness and 
percolation probability can be calculated exactly. The mean cluster 
size is = 2exp[p((i — a)/{l — pa)] — 1 and diverges only at the 
close packing density pcp = l/a. This is confirmed by the behavior of 
the percolation probability. These results should help in judging the 
effectiveness of approximations or simulation methods before they are 
applied to higher dimensions. 

1 Introduction 

One dimensional models have a long history in phase transition studies, go- 
ing back to Ising's solution of the model which bears his name. Such models 
have been found to be free of phase transitions except in some singular cir- 
cumstances (such as zero temperature), but nevertheless keep being studied 
because they are sometimes exactly solvable. Because of this, these models 
can serve as testing grounds, for example, for approximations which are then 
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used in other dimensions. In this spirit, I present here the exact solution of a 
continuum percolation model in one dimension. Like its brethren, this model 
exhibits a phase transition only in some singular circumstances. However, it 
being exactly solvable may yet make it interesting. 

In percolation on a lattice, the one dimensional model is trivial. This 
is not the case for continuum percolation, where the objects connecting to 
each other may occupy arbitrary positions The richness of continuum 
percolation lies in the variety of the binding criterion (which may include 
effects like the shape and spatial distribution of the objects) and the existence 
of interactions. The interplay between the connectivity criterion and the 
interactions make the theory of continuum percolation a very challenging 
field. Therefore, even in one dimension, no completely general results can 
be derived for continuum percolation. However, some specific models can be 
solved. The main such model is the one dimensional version of the extended 
spheres system used some time ago to model microemulsions [0. 

The system consists of + 1 particles on a closed ring of length L. 
The particles interact with each other through a pure hard-core potential 
v{xi, Xj) = v{i,j), i. e., for two particles i and j 



where a is the hard core diameter. 

The connectivity criterion is supplied by the existence of a soft (also called 
permeable) shell of diameter d around each particle. Two particles are bound 
if their shells overlap. Letting p{xi,Xj) = p{i,j) be the probability that two 
particles at Xi and bound, we have 



Naturally, d > a. The percolation transition now arises as a function of 
the density p = {N + 1)/L. The mean cluster size S diverges at a critical 
density pc which signals the sharp (in the thermodynamic limit) appearance 
of an infinite cluster. The order parameter is P{p), the probability that a 
randomly selected particle belongs to this infinite cluster. 

The solution of this model relies on a general mapping I described recently, 
between continuum percolation and a Potts fluid 0. The Potts fluid is a 






X li j I ^ 
3/ 2 j I ^ 
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system of freely moving spins {Xi}^i with s states, which interact with each 
other through a spin- dependent potential V{xi, Xi] xj, Xj), such that 



V{xi,Xi;Xj,Xj) 



U{xi, Xj) if Xi = Xj 
W{xi,Xj) if Xij^Xj . 



[1.3) 



Here U and W are arbitrary functions. 

Furthermore, the spins couple to an external field h{x) through an inter- 
action hamiltonian 



N 



1=1 



^(A) 



:i.4) 
:i.5) 



where 

' s - 1 if A = 1 
-1 if A ^ 1 . 

Every continuum percolation model defined by an interaction v{i,j) and 
a binding criterion p{i,j) can be mapped exactly on a Potts fluid by 
choosing 



U{t,j) = v{i,j) , 
exp[-(3W{i,j)] = q{i,j)exp[-(3v{i,j)] 

where 

(lihj) = '^-p{hj) ■ 
For a constant field, the Potts configuration integral is 



^ = ^ XI / dxi--- dxN exp 



TV 



-^Y.VihJ)+/3hJ2H>^^ 



i=l 



[1.6) 
[1.7) 

[1.8) 



where the sum J2{x^,} is performed over all spin configurations, and /3 = 1/kT 
as usual. 

The magnetization of the Potts fluid is defined as 

1 d\nZ 



M 



and the susceptibihty is 



/3N{s-l) dh 
dM 



X 



dh 
3 



:i.9) 



;i.io) 



The Potts 2-density function is defined as 

INN \ 



\i=i j=i 



According to the general mapping between continuum percolation and 
the Potts fluid, the percolation probability and the mean cluster size are 
obtained as the limits 0, 

P{p) = lim lim lim M , (1.12) 

h^O N^oo s^l 

S{p) = lim lim lini (p < Pc) ■ (1-13) 

Hence the strategy is to calculate M and x for fhs adequate Potts fluid and 
to obtain from them P{p) and S by taking the appropriate limits. Further- 
more, the Potts 2-density function is related to a fundamental quantity in 
percolation, the pair-connectedness function, g\ which is defined as 

g\x^ y) dxdy = Probability of finding two particles in regions 
dx and dy around the positions x and y, such 
that they both belong to the same cluster, 

(1.14) 

where p is the numerical density. This function is related to the mean cluster 
size through the relation (for a translationally invariant system) 

S = l + j^Jdxdyp'g^{x,y) . (1.15) 

The relation to the 2-density function is given by 

g^{x,y) = lim \\p^^\x,p]y,p) - p'^^\x,p-,y, 7])] {p < pc) , (1-16) 

where p,r] ^ 1 and p ^ r], but are otherwise arbitrary. 

For the percolation model defined be Eqs. ( p..l|) and (|1.2| ), the proper 
Potts fiuid is determined by Eq. ( |1.6| ) to be 

exp[-/3f/(x„x,)]^Q(^,j) = I J ^^^Z'^lltl ' (^-^^^ 
exp[-(3W{x^,x,)]^R{^,J) = |° hl^^l^j (1-18) 
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What makes this model solvable is that these functions can take only the 
values and 1 (this condition may be relaxed if the values of d are restricted 
to d < 2a). Mathematically, this system presents several similarities with 
the Takahashi gas , and several parts of the following derivation follow the 
calculations of Takahashi and of Salsburg at al. [§] 



2 Solution of the model 

We assume that the spins are ordered on a closed ring of length L, so that 
the position and the position L are identified. Along this ring, we place 
+ 1 spins, numbered from to A^, so that one of the spins is fixed at the 
position 0. The property which makes such one dimensional models solvable 
is the existence of a canonical ordering of the particles. For a configuration 
{xo, Xi, . . . , a; at}, the canonical ordering consists in labeling the leftmost par- 
ticle (fixed at the position 0), as 0, the one immediately to its right as 1, and 
so on, until the rightmost particle, labeled A^. Thus, 

Q = Xq < Xi < X2 < ■ ■ ■ < xn < L = Xn+I ■ (2.1) 

The definition xn+i = L serves to simplify the notation later on. There are 
exactly (A^+ 1)! configurations which differ from the canonically ordered one 
only by the labels attached to the particles, provided we distinguish clockwise 
numbering from counter-clockwise, i. e., provided the ring is oriented. Since 
all the position variables are integrated upon, each of these configurations 
contributes the same to Z, which is therefore (A^ + 1)! times the contribution 
of the canonically ordered configuration. Hence, 



Z = dxN dxN-i--- I dxi 2^ exp 

Jo Jo Jo rT^i 



{A„} 



TV 

i>j 1=0 



(2.2) 



This property is completely general and is independent of the specific form 
of V{i,j). Let us denote 

h^,,^{\x, - x,\) ^ exp Mr(x„x,)] = I = ^ (2.3) 
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Now, taking into account that xq = xn+i because of the periodic boundary 
conditions, the configuration integral can be written as 



Z = dxN dxN-vJ^ c?xi 2], n A.,A, (|a;i -xjDexp 



N 

(2.4) 

The specifics of the model, Eq. ( |1.18| ), now enter to verify the following 
theorem. 

Theorem (Factorization). For any spin configuration {Am}, 



N 



N-l 



n h^^Ji\^i - = n /A,+i,A.(a;i+i - Xi) . 



(2.5) 



i>j=0 



i=0 



Although simple, this theorem is not entirely trivial because of the pres- 
ence of two scales, a and d, the interplay of which depends on the spin 
configuration {Am}- The theorem is proved in the Appendix. Using this 
result, Eq. ( pl^ ) can now be rewritten as 



N 

x/Ai,Ao(a;i)exp /5/iX!V'(Ai) 

Let us define 



i=l 



(2.6) 



Ga.a,(2/) = /a.a,(2/) exp 1^ [^(A.) + ^(A,)] 



(2.7) 



Then 



Z = [ dxN ■■■ [ dxi V] \Gxo,Xj^iL - xn)G, 

{A„} I 



Ajv,AAr_i(a^Af — XN-i) 

^■■■Gx2M{x2-Xi)Gx,,Xoi^i)^ . (2.8) 
This multiple integral has the form of a Laplace convolution. Denoting 



TA,,A,(a;) = dye--yGx^^xM 



(2.9) 



6 



we have 

roo 



roo 

{Ao,Ai,...,Ajv} 



= Tr(r^+i) . 



(2.10) 



The matrix Txi^Xji'^) now plays the role of an effective transfer matrix for 
this problem. It is an s x s matrix, of the form 



(2.11) 





( ^ 


B 


B ■ 






B 


C 


D ■ 


■ D 




B 


D 


C ■ 


■ D 






D 


D ■ 


■ c) 



where 




A = 


POD 

/ dye 
Jo 


B = 


roo 

/ dye 

Jo 


C = 


roo 

/ dye 
Jo 


D = 


roo 

/ dye 

Jo 



^^M'^^-exp 
to 



-ud+-f3{s - 2)h 



-/3h 



UJ 



Let {Q!i}|=i be the eigenvalues of T. Then 

rdLe-^Z{L)=±ar' 

•'^ i=l 

The eigenvalues are found easily to be 

1 



a2 



A + C+{s-2)D + ^/A 
A + C + {s-2)D - ^/A 



where 



as = ■ ■ . ^ as ^ C - D , 



A = [A-C-{s- 2)Df + 4(s - 1)5^ 



(2.12) 



(2.13) 



(2.14) 
(2.15) 
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Reversing the Laplace transform in Eq. ( p.l3| ) yields 

* 1 ra+ioD 

Z = Y, — :/ du exp[uL + {N + l)\naj{uj)] . 

7=1 27ri J a—ioo 



(2.16) 



In the thermodynamic limit, the integral can be evaluated exactly by steepest 
descent. Let us introduce the notation 



N + 



(2.17) 



Then the maximum of the exponents in Eq. ( p.l6| ) is obtained from the 
condition 



duj 







(2.18) 



which defines the quantity Uj at which 7j(ct;) is maximal. From the usual 
theory of steepest descent 0, we now obtain from Eq. ( p.l6| ) that 



where 



Z = 5:ir,exp [(iV + l)7,(^,)] 



(2.19) 



2n{N + l)^';{uj,) 



1/2 



(2.20) 



The sum in Eq. ( p.l9|) is dominated by 71. To see this, note first that 
since A > 0, we have ai > 02- Furthermore, 



«i — 03 



A-C + sD + ^ 



(2.21) 



Since A = Ce^^^ and h > 0, we have that A — C > 0. Therefore, ai—a^ > 0. 
Hence, for any given u, ai(a') > aj{uj) for j = 2, . . . , s. Since 71 (cj)— 7^(1^7) = 
In (ai/aj), we also have that for any given u, 71(0;) > 'jjiuj) for j = 2, . . . , s. 
Now, at uJi, 7i is maximal, so that 



7i(^i) > 7i(^j) > 7j(^j) 



(2.22) 
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for j = 2, . . . , s. Hence, for — ^ oo 



if, 




1 + E ^ exp [iV7i(^i) - N^A^j)] 

j=2 -"-1 



lim — \nZ 

N^oo N 



7i 



(2.23) 



(2.24) 



From Eq. ( ^ ), we can obtain the maErnetization M, 



M 



d\nZ 



1 



/?(s-l)A^ rf/i 



/5(s- l)iV 
/3{s-l) [dh 



dh 



+ 



' duji 
'dh 



(2.25) 



where we made use of the fact that = (dji/dui)^ by definition of Ui. 

Although Eq. (|2.24| ) is in principle the exact solution of the one dimen- 
sional Potts fluid, it is not explicit enough to be useful. However, we are not 
interested in the Potts fluid itself, but rather in the percolation model. This 
is obtained in the limits s — 1 and h ^ 0, and in these limits all quantities 
can be calculated explicitly. 

To do this, we set s = 1 + e, and calculate all relevant quantities to first 
order only in e (this turns out to be sufficient). Therefore, from Eqs. ( p.l2| ) 
and ( p. 17 ), we find after some algebra that 



7i = uj— + lnA + e 



ph + 



A{A-C + D) 



(2.26) 



where A, B, C, D must be evaluated to zero order in e. From Eq. ( p.l2| ), 
the result for this is 



A 
B 
C 
D 



— exp [— cjcrj 



exp 



-ud — -(3h 



— exp [—ijja — jSh] 

— exp [—ujd — Ph] 



(2.27) 
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Now, 



(3e \ dh 



- — 

^ (3dh 



A{A-C + D) 



(2.28) 



- uJi,e=0 



After some algebra, and with the use of Eq. ( p. 27] ), we end up with 

B ^2 



M = 1 - 



A-C + D 



LJi,e=0 



(2.29) 



After some more algebra, the susceptibility is obtained straightforwardly. 



X 



dM 
'dh 



B^{A + C-D) 
{A-C + Df 



(2.30) 



- aJi,e=0 



In both Eq. (^) and Eq. (g^, Ui is evaluated at e = 0. From Eqs. (^) 
and ( p.27| ), we have that for e = 0, 



7i = ooa — Inuj , 

P 

where p = (A^ + l)/L is the density. Hence from {d'ji/duj)^ 

P 



e=0 



1 — pa 



(2.31) 
0, we obtain 

(2.32) 



The final result is obtained from Eqs. ( p.29[ ) and (|2.30|) by setting e = 0, 
plugging in the value of uJi and taking the limit /i 0. We have 



M P{p) 



1 - hm n 



1) exp 



p{d- 



cr 



1 — pa 



2 exp 



p{d-a) 
1 — pa 



(2.33) 
(2.34) 



If we just set = in the equation for P{p), we obtain P{p) = 0, which 
shows there is no transition. There is however, one exception, obtained at 
p = 1/a = Pep- This is the close packing density in one dimension, and clearly 
we must have a percolating cluster in this case. Indeed, at this density 



lim = exp 



p{d — a) 
1 — pa 



oo 



(2.35) 
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and therefore 



1) exp 



p{d — a) 


r 


1 — pa 





(2.36) 



P=Pcp 



for all values of h; hence, P{pcp) = 1 as expected. Note that for a = (no 
hard core), pcp oo. Eq. ( p^.33| ) becomes in this case 



Pip) = 1 - hm 

h—^0 



1 e 



pd 



(2.37) 



and indeed P{p = oo) = 1. 

The transition is clearer in the expression Eq. (p. 34 ) for the mean cluster 
size. 

p{d — a) 



lim S(p) = lim 2 exp 

P^Pcp P^Pcp 



1 = oo , (2.38) 

1 — pa 

so that indeed the mean cluster size diverges as the density approaches the 
close packing value. 



3 The pair-connectedness function 

Let us now consider the pair connectedness. To this end we need to calculate 
the spin 2-density function, 

INN \ 

P^^Ha;,P;2/,^) = (lZI]5(a;i -a;)5(xj -?/)5A„M^A,,r,) • (3.1) 

\i=l J=i / 

To simplify the algebra let us use the system's translational invariance (one 
can work with Eq. ( |3.lD and obtain translational invariance directly, but this 
adds needlessly to the calculations). We assume that y < x, denote x — y = r 
and set y = and Aq = /i. We now have, therefore, that 

P^^H0,p;r,r7) = p(2)(r;/x,r/) =p^^(5(xfc-r)(5Ao,^5Afefl^ , (3.2) 

where p = {N + 1)/L. From Eq. ( p.8|) , we have that 

p'^ >{r;fi,T]) = 7^y2 dxN--- dxi 
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{Xm} ^ 

X ■■■Gx2Mi^2-xi)Gx^,\o{xi)\ , (3.3) 



where the function Gx^x^ix) is defined in Eq. ( |2.7| ). This is again a Laplace 
convolution and we have 



° ^ fc=l{Ao,Ai,...,Ajv} 

X 5Afc,,?e~'^''/ifc(r; Afc,Ao) , 
where /ifc(r; A^, Aq) is defined by its Laplace transform as 

/•OO . ^ 

/ dre-^^h,ir;X,,Xo)= E Tx,,x,_A^) ■ ■ ■Tx,,Xo{^) = {t') 

•/ r \ \ T 



N 



{Afc-i,...,Ai} 

The sum over the spins yields 

Z ^ 



Afc,Ao 

(3.4) 



^dLe-'^'^-p'^'\r■,^^,v) = J2(T''-'+\u)) e'^^ /.^(r; r/, /i) . (3.5) 
Inverting the Laplace transform in Eqs. ( p.4|) and ( |3.5|) , we obtain finally 



7. , , 1 /-a+ioo 1 ra+ioo ^ ^ 



-lOO 



k=l 

X [r'=(z/) 



(3.6) 



Let us introduce now a basis of eigenvectors of the matrix T{Q, which will 
be denoted {{ui\, (M2I, . . . , {uj \ being associated with the eigenvalue aj 

defined in Eq. ( |2.14| ). Straightforward calculation yields that 

. {B{s — l),ai — A,ai — A, . . . ,ai — A) 

~ [52(s-l)2 + (s-l)(ai-A)2]^/2 ' 
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{B{s - l),a2 - A,a2 - A, . . . ,02 - A) 
[B^{s-iy + {s-l){a2-Ayf^ 

-^(0,1, -1,0,. ..,0) , 
-^(0,1,1, -1,0,. ..,0) , 



(0,l,l,...,-(s-2)) 



(3.7) 



where A, B,ai,a2 are defined in Eqs. (p.l2| ) and ( p.l4| ). Tliese expressions 
can be simplified, however, by setting the magnetic field h = 0, since it plays 
no part in the calculation of the pair connectedness. With this, we have that 
A = C and B = D. From Eqs. ( |2.12|) and (|2.14|) , we now have that 



Hence, 



(«2| 
(«3| 



02 

B{0 



A + {s-l)B 

A-B , 
1 



c 



exp(-C(T) 



^ exp(-Crf) 



1 



1) 



s(s-l] 



1,-1,-1,..., 



^(0,1,-1,0,...,0) 



(3.8) 



1 



s-l){s-2) 



(0,l,l,...,-(s-2)) 



(3.9) 



Finally, we define two vectors of the type (0, 0, . . . , 1, 0, . . . , 0), denoted 
and (r/l such that their single non zero component " 1 " stands at the position 
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/X and r] respectively. With this, we have that 



T 



N-k+l, 



m=l 

s 



(3.10) 



Hence, 



/ du 

J a—ioo 



T 



n=l 



N~k+1. 



uj}\ e 



w(L-r) 



1 ra+ioo ^ 
ZTTt J a—ioo m=l 



E 



a+joo 



where 



[ 2TTi J a 



{fl\Um) {Umlv) 



(3.11) 



(3.12) 



[am(t^)] 

and 7m('-^) = ujL/{N + 1) +lnam(i^) as in Eq. ( |2.17 ). In the hmit N,L ^ oo 
the integral can be evaluated exactly by steepest descent. The calculation is 
identical to the one performed in Eqs. ( p.l6| )- (|2?T^ ). The result is 



ra+ioo 

lim / duj 

N^oo 27li J a—ioo 



T 



N-k+l, 



UJ 



lim Km exp [(A^ + l)'jmiuJm)]Fm,kiuJm, r) 



N^oo 



(3.13) 



m=l 



with Km and uOm defined in Eqs. (|2.20|) and (|2.18|) . As in Eq. (|2.24| ), the term 
7i dominates all the others, so that in the limit oo, we have, exactly. 



ra+ioo 

lim / duj 

J a—ioo 



T 



N-k+l, 



ML-r) 



lim Kiexp[{N + l)-fi{uJi)]Fi^k{uJi,r) 



N^oo 



(3.14) 



From Eq. ( p.24| ), we have that, in the same limit, Z = Ki exp[(A^ + 1)71(071)]. 
Therefore, combining Eqs. (|3.6| ), ( p.lO|) and (|3.14|) , we finally obtain the 
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expression 



p^^\r;fi,r]) = lim p ^ Fi,fc(^i, r)-^ 



N 



N^oo 



k=l 



ra+ioo 






/ du 


T\u) 




' a—ioo 







(/i|Mi)(Mi|r/) 



fc=i 



2Tii 



l+JOO 



n=l 



^{v-ijJi)r 



(3.15) 



In order to calculate the pair connectedness, we need the case ^,1],^ 1 
[see Eq. ( |1.16| )]. For this case, we see directly from the definition of the {ui\, 
Eq. (IJ), that for any A > 2, 



(ni|A) 

(«2|A) 



1 

7^ 



(3.16) 



and for any n > 3, 



(Mn|A) 



(n-2) 



\/(^-2)(n-l) 



if n < A 
if n = A 



(3.17) 



if n> X 



[ J{n-2){n-l) 



Hence, 



^(r;/i,r/) = — / rfz/ ^ (r/|M„) (m„|p) 

■S fc=i ^71"^ Ja-ioo 



ai(c^i) 



(3.18) 

We are interested, however, in the pair connectedness,which involves the 
following difference. 



p(^^(r;p,p) - p(^)(r;p,?7) 



p^ 1 



S 27li Ja 



n=l 



1 k 



,(i'-a;i)r 



(3.19) 
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From Eq. ( p.l6| ), we see that {^\un) — {vl'^n) vanishes identically forn = 1, 2. 
Now, for n > 3, we have that 



ar,\iy 



1 



(3.20) 



and is independent of the value of n. Finally, from Eq. ( p.l7| ), it is easily 
seen that for any fi,!] ^ 1, fi r], we have that 



(3.21) 



n=3 



Therefore, we have that 



P 1 ra+ioo 

S j^—-^ ATTl J a—ioo 



1 k 



(3.22) 



Finally, to obtain g^{r), we take the limit s — 1. In this limit, Eq. (|3.8| ) 
yields 

a;i(co'i) = — exp [— cjicr] , (3.23) 



where, from Eq. ( p.32|) . 



00, = —^ . (3.24) 
1 — pa 

Hence, substituting the values of a„ and ai, we finally obtain , by substituting 
Eq. (|3^) into Eq. ([TTBD , that 



Urn \ [p(^) (r; fi, p) - p^^^ (r; p, r/) 



\ ra+ioo 



P 2vrZ Ja-ioo V'^ 



. (3.25) 



or 



-I oo k / 1.\ 1 i-a+ioo 1 

P f.^lj=0 \3 J Zm J a- loo V 

(3.26) 
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This integral is an immediate inverse Laplace transform. The final result is 
therefore 



1 



oo k 



X Q[r — ka + i{d — a)] exp 



{k-l)\ 

p{r — ka) 
1 — pa 



(3.27) 



where 9(z) is the step function 




(3.28) 



We can understand the meaning of k and j as follows. Given a distance 
r > from the (arbitrary) origin, g\r) is a sum of contributions from sets 
of configurations indexed by k. The derivation shows that in each such 
configuration there are exactly k particles within the interval r. Since the 
potential includes a hard core of diameter a, we expect a condition on k such 
that r > ka. This is indeed contained in the expression Q[r — ka + j{d — cr)], 
since it implies that r > ka + j{d — a) > ka. Next, we note that each set of 
configurations containing k particles in the interval r is further divided into 
subsets indexed by j. The meaning of j is made clear from the condition 
r > ka + j{d — a) = jd + (/c — j)a generated by the step function. This 
implies that j particles out of the k are singlets, i. e., are not connected to 
any other particle in the set. Therefore the step function shows that (7^(r) is 
a sum of separate contributions from all sets of configurations containing k 
particles in the interval r, j of which are not connected to any other particle 
in the set. 

Finally, we can check the expected relation between the mean cluster size 
5* and the pair connectedness. This is worked out most conveniently from 
Eq. (|3.25|) . From this, we have that 



/ dr gHr) 
Jo 



-| oo 



P 



exp 



p{d — a) 
1 — pa 



l/=LxJl 



- 1 



(3.29) 
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From Eq. ( |1.15| ), and remembering that we have calculated g'^{x^y) under 
the assumption that y < x, we now have that 

1 f 

S = l + - j dxdy p'g\x,y) = l + 2p dr g\r) . (3.30) 

p{d — a) 



and therefore 



S = 2 exp 



1 — pa 



(3.31) 



in agreement with the previously obtained result, Eq. ( p.34|) . 



4 Summary 



I have obtained the exact solution of a one-dimensional percolation model 
which includes the effects of interactions among the particles, modeled by a 
hard core of diameter a. The particles are connected to each other if their 
centers are closer than a distance d. The clustering depends on the density p 
of the particles. By mapping this system on a Potts fluid, one can calculate 
the mean cluster size, pair connectedness and percolation probability of this 
system. As with other one dimensional systems, the phase transition oc- 
curs only in some extreme circumstance, in this case, only when the density 
reaches the close-packing value pcp = 1/a. As the density increases towards 
this critical value, the mean cluster size increases as 



S = 2 exp 



p{d — a) 



1 



1 — pa 

and diverges at pcp- Meanwhile, the pair connectedness is given by 



(4.1) 



1 oo fe 

H k=l j=0 



-IV 



p 



[r — ka — j{d — a)] 



k-l 



X 9[r — ka + j{d — a)] exp 



{k-l)\ 

p{r — ka) 
1 — pa 



■ (4.2) 



These results should be particularly useful as a test of approximation 
methods or numerical calculations or simulations. Such methods can be 
checked against the exact results and some idea thus obtained of their reli- 
ability and effectiveness. Once a method has proved to be successful in the 
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one dimensional case, it can be applied to higher dimensions with some hope 
of success. 

From a more general point of view, the present results are a confirmation 
of the power provided by the mapping between continuum percolation and 
the Potts fluid. Recently, Cinlar and Torquato have used renewal theory to 
investigate one dimensional continuum percolation, but without interactions 
{i. e., the ex — *■ limit of the present model) 0. There seems to be no 
simple way of extending this type of direct probability arguments to cover 
interactions as well. It appears that in order to discuss general models, which 
include inter-particle interactions, there is at the moment no method more 
powerful than the mapping with the Potts fluid. 

5 Appendix 

The proof of the factorization theorem is obtained by induction on A^. For 
= 1 the theorem is obvious. Let's assume it is correct for and prove it 
for + 1. Then 

N+l 

n fx.,x,i\xi - 

i>j=0 

N N 
= n h^,xA\^i - ^j\)I[f^N+i,\^ixN+l - Xi) 
i>j=0 i=0 
N-1 N 

= n fx, + l^^{Xi + l - Xi) Yl fxN+l^^{^N+l - Xi) 

i=0 i=0 
N-1 

= n [fx^+iAii^i+l ~ ^i)fxN+i,Xi{^N+l - Xi) fxj^+i,Xj^{xN+l - Xn) , 
i=0 

(A.l) 

where we used the induction assumption to obtain the second equality. 
Let us now consider the various possibilities for every term i. 

1- fxj^+iXi^N+l - Xi) = 1. 

In this case clearly, 

fx,+uxA^i+l ^ Xi)fxj^+i,xXxN+l - Xi) = fx,+-,,xA^i+l - Xi) . (A. 2) 
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2. f\r^+^,\,{^N+i - Xi) = and Xn+i = \i ■ 
Then, from Eq. ( p.3| ), 

fxN+uX^ixN+l - Xi) = Q{xn+i - Xi) = , (A. 3) 

so that, necessarily [see Eq. (|1.1|| , x^+i — Xi < a. However, in the 
canonical ordering, Xi < Xj+i < • • ■ < a;Ar, so that 

xat+i — xn < xn+1 ~ Xi < a < d . (A. 4) 

Hence, 

/Ajv+i,Ajv(a^Af+i — a^Af) = , (A. 5) 

even if Aat+i Xn, because a < d. Therefore, 

= f\i+i,Xi{^i+l ^ ^dfxN+iXi^N+l — Xi)f\j^^^^\j^{xN+l — Xn) 
= f\i+i,Xi{Xi+l-Xi)fxj^+i,Xj^{xN+i-XN) ■ (A. 6) 

3. fx,^+uxX^N+i - Xi) = and Aat+i ^ Xi . 
Then [see Eq. 

fxN+l^^{xN+l - Xi) = R{xn+i - Xi) = . (A. 7) 

and, necessarily [see Eq. ([lT[| , xn+i — Xi < d. Then there must exist 
some j,i<j<N + l, such that 

Aat+i = Afe and Aat+i ^ Xj for all j < k < +N + 1 , (A. 8) 

i. e., j is the closest spin to Aat+i which differs from it. In particular, 
Atv+i = Aj+i, so that 

fx,+^,x,ixj+i - Xj) = R{xj+i - Xj) . (A.9) 

Because of the canonical ordering, 

Xj+i — Xj < Xj+i — Xi < xn+1 — Xi < d . (A. 10) 

Hence 

R{xj+i-Xj)=0 . (A.ll) 

Therefore, 

= fx,+i,xA^i+l ~ ^i)fxN+i,xA^N+l - Xi)fx^+^,Xj{Xj+l - Xj) 

= /A,+iA(a;m -a;i)/A,+i,A,(a;j+i -Xj) . (A.12) 
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From Eqs. ( [A.2| )-( [702D , it follows that fxr^^j^^xX^N+i — Xi) makes no dif- 
ference anywhere in the product Jlill^ /A,;+i,A,(a;i+i - Xi)fxf,+^x{^N+i - Xi). 
It can therefore be dropped out of every term i. Hence, 



n [f\+l^^{Xi+l - Xi)f\^^^,\,{xN+l - Xi)\ fx,,+„X^{xN+l - Xn) 



i=l 



N~l 

n [fx.+uxX^i+i - Xi)] fx^+„x^{xN+i - Xn) 

i=l 

N 



Ylfx^+ixi^i+i - Xi) 



(A.13) 



i=l 



Referring back to Eq.( [A.l[ ), we see that this proves the theorem. 
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